Static friction on the fly: velocity depinning transitions of lubricants in motion 
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The dragging velocity of a model solid lubricant confined between sliding periodic substrates 
exhibits a phase transition between two regimes, respectively with quantized and with continuous 
lubricant center-of-mass velocity. The transition, occurring for increasing external driving force F ex t 
acting on the lubricant, displays a large hysteresis, and has the features of depinning transitions in 
static friction, only taking place on the fly. Although different in nature, this phenomenon appears 
isomorphic to a static Aubry depinning transition in a Frenkel-Kontorova model, the role of particles 
now taken by the moving kinks of the lubricant-substrate interface. We suggest a possible realization 
in 2D optical lattice experiments. 

PACS numbers: 68.35.Af, 05.45.Yv, 62.25.+g, 62.20.Qp, 81.40.Pq, 46.55.+d 



One of the current areas of development in the modern 
understanding of sliding friction is the depinning tran- 
sition between substrate and slider in static friction 
A slider initially pinned on a substrate requires a finite 
force in order to move, whereas an unpinned one will 
move even under an infinitesimal force, ft is generally 
known that ideal crystalline sliders commensurate with 
periodic substrates are always pinned, whereas incom- 
mensurate ones exhibit as a function of material param- 
eters a sharp transition between a rigid, pinned state for 
soft sliders on strongly attractive substrates and a freely 
sliding state for hard sliders on weakly attractive per- 
fect substrates. The standard paradigm for that transi- 
tion is the Aubry transition in the time-honored one- 
dimensional (ID) Frenkel-Kontorova (FK) model 0, 0|- 
Once sliders are set into motion, however, one would gen- 
erally expect no further rigidity or pinning to persist or 
to appear, ft was thus rather surprising to find that 
simple systems like the 3-lengthscale driven ID model of 
Ref. [B| , inspired by the tribological problem of two slid- 
ing crystal surfaces with a thin solid lubricant layer in 
between, exhibits a novel kind of rigidity, now around a 
state of dynamical motion, where the lubricant center- 
of-mass (CM) velocity is robustly locked onto quantized 
plateau values 0, 0> Hi ■ The dynamical pinning of the 
lubricant's motion onto a rigidly quantized relative ve- 
locity state has been understood in terms of kinks of the 
lubricant being set into motion by the shear due to the 
moving surfaces. This kind of kink dragging can be ar- 
gued to represent a rather general mechanism, possibly 
at play also in more realistic two-dimensional situations. 

One such situation could be realized e.g., by attempt- 
ing to slide a rare gas physisorbed layer on a crystalline 
surface, by means of an external dragging agent, such as 
an AFM tip, or an optical tweezer, or a 2D optical lattice. 
In addition, the substrate could be oscillated in a Quartz 



Crystal Microbalance (QCM), the inertial force acting as 
an (infinite wavelength) dragging agent. In either case, 
the external potential would in reality drag the solitons 
or discommensurations formed by the adsorbate with the 
substrate - it would drag the 2D Moire patterns. While 
such a class of phenomena has not yet been explored, it 
is potentially quite interesting, and its extent and con- 
sequences deserve a full theoretical understanding. In 
this Letter we present a close analysis revealing that the 
lubricant velocity-pinned state shares many more of the 
characters typical of ordinary pinning in static friction, 
including hysteresis against a depinning external force 
-Fext, and a genuine depinning transition as a function of 
lubricant stiffness, that appears a dynamical isomorph to 
the Aubry transition of static friction. 

Starting with a chain of harmonically interacting par- 
ticles Q (the lubricant layer), the Hamiltonian f 
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contains the interaction with the two mutually sliding 
substrates of sinusoidal amplitudes C/+ and U— and pe- 
riods a+ and a_, different from the rest length ao of 
the springs, which sets a third length scale. Introduc- 
ing dimensionless length ratios r± = a±/ao we assume 
r_ > max(r + , r7 ), whereby the (+) slider has the clos- 
est periodicity to the lubricant, the (— ) slider the fur- 
thest. To study depinning, we here apply a uniform 
external force F ext to all chain particles. The infinite 
chain size is managed - in the general incommensurate 
case - by means of periodic boundary conditions (PBC) 
and finite-size scaling 0, Sliding the substrates at 
constant relative velocity w cx t, the lubricant equations 



of motion are integrated by a standard Runge-Kutta al- 
gorithm, including a phenomenological viscous friction 



force —27 (±i 



t) = -27 {Xi 



,) accounting for 



various sources of dissipation. As shown earlier [(J 0], 
for F cxt = and a wide range of parameters, the sliding 
chain reaches a dynamical stationary state characterized 
by a quantized plateau velocity 



^plateau /^ext — 1 
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depending solely on the incommensurability ratio r+ , in- 
dependent of r_, K, 7, w cx t, and even of U-/U+. In this 
quantized state of motion the array of topological solitons 
(kinks for r + > 1 or antikinks for r+ < 1) formed by the 
chain with the least mismatched (+) substrate is rigidly 
dragged at velocity v G xt by the most mismatched (— ) 
substrate. This dynamical pinning of the kink velocity 
to the moving (— ) slider suggests a more complete anal- 
ogy to the static pinned state of particles in the standard 
FK model, where a finite force F ext > F FK is required 
to slide the chain over a periodic substrate. Depend- 
ing on the ratio 9 FK between the substrate and chain 
periodicities, static pinning in the FK model occurs for 
arbitrary K for rational 9 FK , whereas irrational 9 FK, s 
show a depinning (Aubry) transition taking place at a 
critical chain stiffness i^Aubry, beyond which free sliding 
is induced by arbitrarily small F ext [H, 0] • In the present 
slider-lubricant-slider problem, the kinks formed by the 
lubricant over the (+) slider form a ID lattice with an 
average distance = a + /(r + — 1), implying a kink cov- 
erage 9 = a^/dk = r_ (1 — r^ 1 ) of the (— ) substrate (for 
antikinks, 9 < 0). If the kinks here play the role of the 
particles in the FK model, we should be able, depend- 
ing on the relative commensurability 9, to observe either 
indefinite pinning (i.e., quantized lubricant velocity) in 
the commensurate case, or else an Aubry-like depinning 
transition in the incommensurate case. We find that this 
is what happens, moreover with a clearly outlined iso- 
morphism of the dynamical kink pinning problem to the 
static particle pinning of the FK model. 

In order to highlight this deep similarity with the static 
case, beginning from a quantized-velocity state, we start 
off by studying the motion of the chain under the action 
of an adiabatically cycled external force F eX f We select 
two prototypical kink coverages, one where the kink lat- 
tice is commensurate with the (— ) dragging slider (6=1, 
realized by r + = <fi, r_ = <fi 2 , with </> = (1 + y/o)/2, the 
golden mean), and another incommensurate (9 = 4>, ob- 
tained with r+ = <f), r- = (f) 3 ). The exactly quantized 
CM velocity plateau at F ex t = suggests null response 
to perturbations, so that any < F cxt < F^" should 
have no effect whatsoever on u cm ( "incompressibility" ) . 

Figure [T] displays the result of cycling the external 
force F ext up and down in small steps, for 9 = 1, 
7 (K m)^ 1 / 2 -C 1 (underdamped regime), two different 
i> ex t, and a chain stiffness K within the velocity plateau. 
It displays a clear incompressibility and a hysteretic loop, 
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FIG. 1: (Color online) Hysteresis in the v cln — F ex t charac- 
teristics for (r+,r_) = ((f), (f> 2 ). Throughout, our units are 
m = 1, a+ = 1, F+ — 2irU+/a+ = 1, K is measured in 
F+/a+, 7 in J mF+ /a+ , velocities in ^/a+F+/m, and we 
choose F- = 2-kU-/ci- — F+. Q The behavior is shown for 
fast (t) ex t = 0.1, upper panel) and slow (f cx t = 0.01, lower 
panel) drive. Adiabatic increase and decrease of F GXt is de- 
noted by triangles and circles, respectively. A characteristic 
multi-step feature appears when decreasing adiabatically F ex t- 
Here 7 = 0.1 and K = 4. 



analogous to the depinning transition in static friction [l[ , 
and in the FK model 0]. For increasing F GX t the lubri- 
cant average velocity w cm is discontinuous at F c + 1 , with a 
dynamical depinning, a finite jump Av and a bistable be- 
havior. When F oxt is decreased back, the depinned state 
survives down to F^ < F+l, where quantized sliding is 
retrieved, as shown in Fig.[TJ The hysteretic behavior ob- 
served under the application of a negative -Fcxt is similar 
to that for F oxt > 0, with corresponding < < 0. 
Viewing the kinks of our model as the particles of a stan- 
dard FK model, the 9 = 1 case is isomorphic to the fully 
commensurate 9 FK = 1 FK static case, where pinning 
holds for all K' s. Confirming that, when we counterbal- 
ance the average frictional force F w = —27 (u p i a teau — ^w) 
with an equal and opposite external force F cxt = — F w 
such that no net force acts on the kink lattice, we find 
that the plateau extends out to arbitrarily large K, see 
Fig. [D(a). The observed phenomenology is quite generic: 
in suitable parameters ranges we find similar hysteretic 
plateaus for all values of r + and r_ investigated. More- 
over, different choices of r_ = 9/(1 — r^ 1 ), with rational 
9^1 lead to plateaus which also extend to arbitrarily 
large K . A generic rational 9 > 1 implies defects in the 
perfectly commensurate (9 = 1) lattice of kinks, in the 
form of kinks of the kink lattice — hierarchical excita- 
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FIG. 2: (Color online) (a) The critical depinning (F+ T , F~ l ) 
and re-pinning (F^, F^ ) forces, as functions of the chain 
stiffness K, for (r+,r_) = (4>,4> 2 ) (i.e. commensurate 9 = 1, 
solid symbols, pinning extending to infinitely large K) and for 
(r+,r_) = ((f), (f) 3 ) (i.e. incommensurate 8 = <^>, open symbols) . 
In all calculations v cx t = 0.1, 7 = 0.1. For 9 = 1: (b) the 
width F c +T — F^~* of the bistability region as a function of 
(K, — K) (crosses), compared to a fitted power- law F+ T - 
F+ l = B {K,-K) a , with a ~ 2.7, K* ~ 135; (c) the large- JiT 
depinning force, compared to the expected power-law (F c + T + 
F w ) = Aif" 1 , Eq. ©. 



tions, fcinfe 0/ fcinfe — which, being commensurate with 
(— ) substrate, resist depinning for arbitrarily large K. 

On the contrary, we expect — and find, see Fig. [2ja), 
for 9 = <j> — that irrational 9 should be associated 
with a genuine Aubry transition on the fly of the in- 
commensurate lattice of kinks, with the plateau disap- 
pearing at ^Aubry' The transition takes place at a finite 

K = if A ut,ry > wnere the width F+T — of the velocity 
plateau as a function of Fe X t shrinks to zero. Of course, 
we expect to be a complicated function of the 

original model parameters, generally not coincident with 
the value characteristic of the FK model for 9 FK = 9. 

The pinning force (F^ + F w ) decreases like K^ 1 
for large K, see Fig. EJc). In this limit, the kink lat- 
tice becomes increasingly faint, since the correspond- 
ing particle-density modulation amplitude drops as K . 
Quantitatively, the small displacements of the particles 
away from ideal lattice positions Xi — tag due to inter- 
action with the (+) potential are described in terms of 
a hull function, whose explicit expression is known [Iol | 
for large K > F + /a + (with F + = 2nU + /a + ). These 
displaced positions, once substituted into the U— term 
of Eq. ([JJ), yield a sinusoidal oscillation as a function 
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FIG. 3: (Color online) Time evolution of the positions Xi — 
Wpiateau-t of 7 contiguous chain particles, t is measured in units 
of y/ ma+ 1 F + . The plots refer to three dynamical regimes 
observed along the adiabatic decrease of F ex t, at the points 
indicated by (3a)-(3c) in Fig. Q] (a) free sliding at F oxt = 
0.07, (b) dynamic stick-slip at F ext = 0.045, and (c) quantized 
sliding state at F ext = 0.02. Same parameters as in Fig. [T] 
with Wext =0.1. 



of the translation v cx tt of the (— ) substrate relative to 
the kink lattice. This energy oscillation corresponds to 
a force whose maximum amplitude equals the minimum 
force that must be applied to the chain to have the kinks 
overcome the barrier and initiate sliding. This critical 
force per particle amounts to 
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and is drawn in Fig. [2fc) for comparison. The ~ 20% 
discrepancy with the observed F+ is due to the neglect 
of the displacements induced by the (— ) potential, which 
enhance the kinks amplitude in this 9 = 1 geometry. 

The self-pinning at the (— ) minima is also at the root 
of the bistability found at moderate K: self-enhanced 
trapped kinks resist to a large F ox t up to F+* , but 
as soon as kinks wash away in the large- F ext unpinned 
sliding state, F oxt must be reduced to smaller be- 
fore kinks reconstruct and then re-pin. The transi- 
tion boundaries F^ , F c + ^, and At; are nontrivial func- 
tions of the parameters. In particular, as shown in 
Fig. Ufa), for large K the hysteretic behavior extends 
up to K = ~ 135, see power-law fit in Fig. [2]^b). 
For K > K*, F+ T = F+ l = F+ , and the depinning 
transition is continuous, characterized by a critical be- 
havior (v cm - Wpiatoau) oc (F GXt - F+) 1 / 2 . On approach- 
ing the critical point K*, the jump Av characterizing the 
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FIG. 4: (Color online) Hysteresis, as u cxt is cycled, in the 
underdamped regime. 



first-order depinning transition vanishes, and a continu- 
ous transition line originates there. In the moderate- if 
hysteretic regime we find that, in a wide range of param- 
eters, depinning occurs through a discontinuous jump to 
a quasi- free sliding regime F ext > -F+l, see Fig. [IJa), 
characterized by aperiodic single-particle motion, super- 
posed to a drift approaching the free-particle limit veloc- 
ity F cxt /(2j). The reverse route, upon decreasing F ext 
to F^, yields some insight in the re-pinning mechanism. 
Depending on the values of K and u ex t, dynamical re- 
pinning to the plateau-state takes place through a vari- 
ety of mechanisms, from intermittencies with a well de- 
fined periodicity, to more chaotic and irregular jumps. 
Despite that variety, for intermediate K we consistently 
observe one or more multi-step sliding regimes as F ext 
decreases within the [-F+- 1 , F+T] hysteresis window, con- 
sistently with the results known for the sine-Gordon ll| , 
and FK [3] models. Our intermediate- i 7 ^ regime, illus- 
trated in Fig.[3jb), is reminiscent of intermittencies rep- 
resenting a dynamic stick-slip occurring on a local scale, 
clearly seen by plotting the particle trajectories in the 
reference frame which slides at u p i a tcau- In fact, detailed 
global analysis of this intermediate- F ext dynamics, based 
on simulations with r + = 1 + e (with e <C 1) where kinks 
are well separated, shows that the kinks pin to the (— ) 
lattice along most of the chain, with a few kink-antikink 
defects travelling along the chain. These defective sec- 
tions move along the chain at a characteristic speed de- 
pending on the model parameters. The nature of these 
de-pinning mobile defects is reminiscent of those found 
in the ordinary FK model Q . Depending on K and F cxt , 
single- or multiple-kink defects arise: transitions between 
intermediate phases characterized by different numbers of 
travelling-kink defects are visible in Fig. [1] 

In practice, the experimental realization of an equal 
force acting on each particle could be provided, for in- 
stance, by the inertial forces in a QCM experiment. In 
addition, the hysteretic depinning and re-pinning occurs 
even when different tunable parameters are varied adia- 
batically, for instance the driving velocity iw, see Fig. SI 
as long as we are in the underdamped limit [131 ]. Thus, in 
a concrete lab experiment, cycling quantities such as u ex t 



should lead to leaving/recovering the plateau state, simi- 
lar to cycling F cxt . As for the dragging agent, the typical 
force that a laser provides in an optical lattice of alkali 
atoms is F m 10~ 9 pN [14]. This is very small, but we 
checked that our phenomenology still survives even when 
F— is orders of magnitude smaller than F + . Larger (in- 
ertial) forces, mu> 2 a ~ 10 _5 pN per adsorbate atom, are 
accessible in a standard QCM oscillating at uj/2tt ~ 15 
MHz with amplitudes a ~ 100 A. The sinusoidal varia- 
tion of the dragging velocity occurs here between and 
±lm/s on a much longer time-scale than the typical ad- 
sorbate vibrational frequencies (~ 1 THz). The inertial 
dragging force couples to the local density fluctuations, 
and should drag the solitons exactly like a commensurate 
(— ) slider would do. 

As for thermal effects, they favor unpinning, but as 
long as ksT is much smaller than the cost A for the 
creation of the relevant depinning defects, the physics 
does not change much, as shown by numerical simulations 
of similar models. [TBI] 

Our overall picture confirms a striking similarity of the 
dynamical kink-depinning transition of a lubricant under 
shear to the usual depinning in static friction, and in the 
static FK model, to which the large- if theory of Eq. ([3|) 
establishes a quantitative correspondence. These findings 
should also be amenable to experimental verification. 
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